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Abstract 

The idea that the pomeron has partonic structure similar to any other hadron 
has been given strong support by recent measurements of the diffractive structure 
function at HERA. We present a detailed theoretical analysis of the diffractive 
structure function under the assumption that the diffractive cross section can be 
factorized into a pomeron emission factor and the deep inelastic scattering cross 
section of the pomeron. We pay particular attention to the kinematic correlations 
implied by this picture, and suggest the measurement of an angular correlation 
which should provide a first test of the whole picture. We also present two simple 
phenomenological models for the quark and gluon structure of the pomeron, which 
are consistent with various theoretical ideas and which give equally good fits to 
recent measurements by the HI collaboration, when combined with the pomeron 
emission factor of Donnachie and Landshoff. We predict that a large fraction of 
diffractive deep inelastic events will contain charm, and discuss how improved data 
will be able to distinguish the models. 



1 Introduction 



Recent measurements at HERA have indicated that a significant fraction of deep inelastic 
electron-proton scattering events have a final state with a large rapidity gap between 
the proton beam direction and the observed final state particles H|. The lack of any 
hadronic activity around the proton beam direction and the mismatch between the initial- 
state and observed final-state energy requires the proton (deflected only by a small angle 
and therefore outside the rapidity coverage of the detectors) to be in the final state, still 
carrying a large fraction of its incident momentum. These events with a remnant proton 
in the final state are classified as diffractive scattering (DS). In analogy to the conventional 
deep inelastic scattering (DIS) cross section, the differential cross section for DS can be 
written as 



da DS 4yra 2 



2 

da p dtdxdQ 2 xQ 4 I 1 ~ V + 2[1 + R^ix, Q 2 , a p , t)} J Q ^ ^ t} ' (1) 

The final state configuration of these events suggests that they are caused by a deep 
inelastic scattering of an uncharged and uncoloured object, which was emitted from the 
proton beforehand, Fig. |l|. From the kinematical distribution of the diffractive events it 
seems most likely that this object is the pomeron, which so far has only been observed 
and identified by its t-channel trajectory || in the full hadron-hadron cross section^] The 
idea that the pomeron has hard partonic constituents was first proposed by Ingelman 
and Schlein ||, and given strong support by the hadron collider experiments of the UA8 
collaboration || [7|. 

If this interpretation is correct, then one would expect that the diffractive cross section 
could be factorized into a piece corresponding to the emission of an uncharged, colourless 
pomeron from the proton and a piece corresponding to a hard scattering off the partonic 
constituents of the pomeron: 



da DS 47ra 2 



VP + oh , Jf ^ Tvj | F 2 p (z,Q 2 ,t)f(a p ,t), (2) 



da p dtdzdQ 2 zQ A \ ^ 2[1 + R p (z,Q 2 ,t)] 

where z = z(x,Q 2 ,a p ,t) is the fraction of pomeron momentum carried by the struck 
parton and 

yp = — (3) 

is the virtual photon energy 'seen' by the pomeron. F^iz, Q 2 , t) denotes the DIS structure 
function of the pomeron and f(a p ,t) the probability for the emission of a pomeron with 
momentum fraction a p and t-channel momentum t off a proton. 

1 There has been recent evidence j| for a glueball candidate at M = 1.9 GeV, which lies on the timelike 
continuation of the pomeron trajectory and has the correct quantum numbers predicted by Regge theory. 
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A common but only approximately correct way of parametrizing this factorization 
property is to write the diffractive structure function as the product of an emission factor 
and the deep inelastic structure function of the pomeron 0. [| : 

F 2 DS (x,Q 2 ,a p ,t) = F 2 p (z,Q 2 ,t) f(a p ,t). (4) 

Note that the R and R p functions cannot be related in such a simple manner. We 
will discuss various tests of the factorizability of the cross section and investigate the 
applicability of the factorization at the level of structure functions (Q). 

Due to the lack of experimental information on the remnant proton, a complete kine- 
matical reconstruction of diffractive scattering events is not possible at present. Both 
parameters describing the pomeron emission (a p and t) can only estimated indirectly or 
have to be integrated out. The kinematical parameter z, describing the fraction of the 
pomeron's light-cone momentum 'seen' by the virtual photon, can be (up to a small un- 
certainty) obtained by measuring the invariant mass of the hadronic system X in Fig. [I]. 
Since to a good approximation a p = x/z, the a p dependence of F 2 DS can be inferred and 
compared with theoretical predictions for the function /. In this approxmiation, the fac- 
torizability of the structure function (|4]) becomes an exact statement following from the 
factorizability of the cross section (|^). We will discuss the reconstruction of the kinematics 
and the uncertainty on a p caused by the lack of kinematical information on the remnant 
proton in Section 0. 

If the object struck by the virtual photon in diffractive deep inelastic scattering is 
indeed the same pomeron which controls the high energy behaviour of hadronic scattering 
amplitudes, then its basic properties and in particular its coupling to the proton are 
already known. For example, Donnachie and Landshoff give a simple form for f(a p ,t) || 
which they derive from an essentially nonperturbativeQ model \TU\ of pomeron exchange 
dynamics in terms of Regge amplitudes: 

q/,2 

f(a P ,t) = — 2 [F 1 (t)} 2 al- 2a ^. (5) 

The Dirac form factor of the proton entering in f(a p , t) is well known from low-energy ep 
scattering fT3|| : 

4M 2 - 2.8t / 1 



Fl{t) ~ AM 2 — t \l-t/(0.7GeV 2 " ' ' (6) 




2 There have been several recent attempts to derive a perturbative formulation of the pomeron. These 
approaches fill , all based on the BFKL equation will not be discussed in the context of this paper, 
as there is insufficient conclusive evidence at present for the applicability of the BFKL equation in the 
kinematic range covered at HERA. In the following discussion, we will always assume f(a p , t) (as for any 
other hadron-hadron interaction at low invariant momentum transfer) to represent a nonperturbative 
coupling of pomerons to the proton, which can be determined from the experimental data. 
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where M is the proton mass, whereas the pomeron coupling strength to quarks^ b ~ 
1.8 GeV -1 and the pomeron trajectory 

a(t) = l + e + a't with e = 0.086, a' = 0.25 GeV" 2 (7) 

are tuned to explain a wide range of experimental results in elastic pp, pp, and up scat- 
tering [Q]. Other similar forms for / have been proposed in the literature, see for example 
Ref. ||, but the differences are not crucial to the present discussion. 

The above picture has recently been given strong support by a detailed analysis of 
diffractive deep inelastic scattering events by the HI collaboration at HERA p^| . Their 
principle conclusions are: (i) the Q 2 dependence of F® s is consistent with scattering off 
pointlike objects, (ii) the factorization of the diffractive structure function into pieces 
which depend separately on z and a p , Eq. (||), is observed, (iii) the a p dependence of 
/ is consistent with that predicted by Donnachie and Landshoff, i.e. ~ a^" 201 ^, and 
(iv) the pomeron structure function F 2 P is 'hard', i.e. the pointlike constituents carry 
a significant fraction of the pomeron's momentum on average. Not yet determined ex- 
perimentally are: (i) the 'nature' of these hard constituents (i.e. whether the pomeron 
predominantly consists of quarks or of gluons), (ii) the explicit t-dependence of F^ 3 pre- 
dicted by Eqs. (|]J6|J7|), (hi) the pomeron flux factor (i.e. a possible overall normalization 
factor in Eq. (|5|)), (iv) the kinematical distribution of the remnant protons, and (v) the 
magnitude of the i?-factor (its impact on the HI results has been shown to be less than 
17% 0). 

It is these latter issues that are the subject of the present study. Having established 
experimentally that the overall picture is consistent with Fig. [I], the next task is to ask 
more detailed questions. We have already mentioned that the kinematic variable a p 
cannot at present be measured directly, and so it is important to study the uncertainty 
which is introduced when it is reconstructed from observed quantities. The kinematical 
constraints implicit in Fig. [I] also lead to small but non-negligible correlations between 
the final state electron and proton. The magnitude of these correlations depends on the 
form of / and _F 2 P . Finally we shall investigate the z and Q 2 dependence of _F 2 P itself. The 
HI data already contain a significant amount of information. 

In particular, we shall show that the data can be understood either in terms of a 
model in which the pomeron predominantly consists of gluons at a scale Ql = 2 GeV 2 , 
with a small admixture of quarks such that a momentum sum rule is satisfied, provided 
that the pomeron flux factor is rescaled, or in terms of a model in which the pomeron, 
like the photon, has pointlike couplings to quarks such that the structure function has an 
additional 'direct' contribution. In the case of the latter model, the pomeron flux factor 

3 We use the notation b rather than (3 to avoid confusion with the kinematic variable introduced in 
Section ||. 
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requires no further rescaling. We shall discuss in the some detail the different z and Q 2 
behaviours expected in the two approaches. A common feature is the prediction of a 
sizeable charm quark contribution to diffractive deep inelastic scattering. 

The paper is organized as follows. In the following section we study the kinematics of 
diffractive deep inelastic scattering, as implied by Fig. |], in some detail. In Section |3| we 
present predictions for a particular kinematic correlation which should be straightforward 
to measure and which will provide a stringent test of the pomeron picture. In Section |] we 
discuss models for the parton structure of the pomeron, and the corresponding predictions 
for the z and Q 2 dependence of the pomeron structure function F 2 P . Our predictions are 
compared with the experimental data from HI in Section ^. Finally, Section |6| contains 
our conclusions. 



2 Kinematics of electron-pomeron deep inelastic scat- 
tering 

2.1 Reconstruction of the kinematical invariants 

To reconstruct all kinematical parameters in (fj), it is sufficient to measure the momenta 
of the outgoing electron (52) and the remnant proton (p'). It is convenient to parametrize 
these in a Sudakov decomposition using two lightlike vectors directed along the beam and 
a spacelike transverse vector. Since we are ignoring the electron mass we can use the 
incoming electron momentum qi for one of the lightlike vectors. For the other, we define 

M 2 

where s = (p + qi) 2 , p 2 = M 2 and, by construction, p 2 = 0. Hence we can write 

q 2 = Ap + Bqi + q T 

p > = Cp + D qi + k T , (9) 

which implies 

q = -Ap + (1 - B)q 1 - q T 

( M 2 \ 

k = (l-C)p+[^- ] ^-Djq 1 -k T . (10) 

The eight degrees of freedom in (^J) are reduced to five by requiring that q\ = 0, p' 2 = M 2 
and disregarding an overall azimuthal angle. The next step is to relate these to more 
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familiar deep inelastic and diffract ive variables. The electron is described by the usual 
two DIS variables x and Q 2 , and three additional parameters define the proton: 



a. 



fraction of longitudinal momentum transferred to the pomeron, 
t-channel invariant momentum transfer to the pomeron, 
angle between the outgoing electron and outgoing proton 
in the plane transverse to the beam direction. 



In terms of Lorentz invariants, 



Q 2 



-q , x 



Q 2 



t — k 



a. 



k ■ qi 



(11) 
(12) 

(13) 



(14) 



2p ■ q ' P ■ qi 

Some straightforward algebra then gives the result for the photon and pomeron momenta: 

' Q 2 \ 



k 



s-M 2 ) 

t - a n M 2 



a p p + 



s-M 2 



where 



Q 2 l- 



p + 

qi - k T 

_Q 2 

x(s 



1 M 2 
- + - 

X s 



M 2 



Qi 



- Qt 



(15) 



M 2 Q 2 



COS ( 



'ep 



-t(l- 

qr ■ k 2 



OL, 



-M 2 ) 
oclM 2 



M 2 f 



(16) 



As already mentioned, neither a p or t are directly measured. An additional constraint 
can however be obtained by measuring the mass of the final state in the r y*{q)P{k) — > X 
hard scattering, i.e. M\ = (q + k) 2 . In analogy with the usual Bjorken x variable, 
Eq. (|i~4"l), we introduce 

Q 2 Q 2 (m 

Z 2q-k M 2 x + Q 2 -f 1 } 

Substituting the expressions for q and k from Eq. ( ]T5] ) then gives the required constraint: 

1 a p 2a p M 2 - 1 
z x s — M 2 



2 cost; 



'ep 



{-t(l -Op) -a 2 p M 2 } Jl 



Q 2 
x(s - M 2 ) 



M 2 Q 2 



M 2 f 
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From Eq. (||) we see that large values of \t\ are expected to be heavily suppressed, 
and this is consistent with the fact that no final-state protons are observed outside the 
beam pipe. It is therefore a reasonable first approximation to set t = in the kinematical 
relations. With t = M 2 = 0, Eq. (|l8f) becomes 

i = 2E =► a p = - 1 (19) 

Z X z 

with the interpretation that the momentum fraction of the quark in the proton (x) is 
simply the product of the momentum fraction of the quark in the pomeron (z) and the 
momentum fraction of the pomeron in the proton (a p ). Note that in this approximation 

z = Q 2 /(Q 2 + M 2 x ). (20) 

In this way, the parameter a p is easily determined from measured quantities. 

It is important to stress, however, that the corrections to fllQ]) are not obviously negli- 
gible. In particular, we note that the terms of order \f—tjQ and M/Q may not be small, 
while corrections to ( p0|) start at order tj{M\ + Q 2 ) and will therefore be ignored in the 
following. In practice, the necessity to have a large rapidity gap in order to distinguish the 
diffractive events requires the pomeron to be slowly moving in the laboratory frame, and 
consequently a p <l. In this limit, including the most important subleading corrections 
gives 

x 



CXp 



xs Q 



l + 2yl — cos(J) ep 



(21) 



This result shows that the distribution in the angle (f) ep will not be uniform in general. 
For any non-zero t, and at fixed z, x and Q 2 , a p varies with <f> ep . Since the diffractive 
structure function is a steeply falling function of a p , Eqs. (fUD, the effect can be quite 
large. This effect will be studied in greater detail below, and in particular the implications 
for angular correlations between the outgoing electron and the remnant proton will be 
elaborated in Section ||. 



2.2 Estimates for the systematic uncertainties at HERA 

The dependence of a p on the presently unmeasurable angle 4> ep gives rise to a system- 
atic uncertainty on reconstructing the variables yp and a p which appear in (0). In this 
section we attempt to quantify these uncertainties in order to test the validity of the 
appr oximat ions 

ot P « -, yp&y (22) 
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used to extract F<[{z, Q 2 ) from the HERA data 11510. We will also test the factorizability 
of the diffractive structure function (f|). 

For any parametrization of f(a p , t) which has a similar a p dependence to (|5|), one ob- 
tains a diffractive cross section which decreases steeply with a p . Therefore the correction 
will not average out over all angles <f) ep , rather these corrections will accumulate to 
give a non-zero average deviation from (|22l) . The relative deviation of a p from x/z is 
given by 

a p — x/z n / 1 Q 2 Z\f—i 



x/z V xs Q 

while the relative deviation of yp from y has a similar form: 



2\ 1 — — cos cj) ep , (23) 



y P -y x Q 2 z^i a p -x/z 
= 1 ^ -2W 1 -— cos ep = 

y za p V xs Q x/z 



(24) 



Since (p ep and t are not directly measured, we define the expectation value of the deviation 
to be the weighted average over all angles and values of t|: 



r 2n 



•a x/z\ J dt Jo ^ P (^j^)f(a p ,t) 
a ?- x/z \(x,z,Q 2 ) = ^4 V X/Z I , (25) 



x/z 



— oo 



dt / d(f) ep f(a p ,t) 



Vp ~ V \ (x,z,Q 2 ) = -l^±\ (x ^) (26) 



\ y I \ x/z 

which becomes, for any f(a p ,t) with a similar form to 



' a p -x/z\ 2 / Q 2 \ ^ 2 



dt (-f)(l-2a(f))/(^,f 



) = 2 1 - f- ^ / . (27) 



x/z / ' ' \ xs I Q 



Figure ^|(a) shows this systematic deviation for the DL parametrization of f(a p , t) (|5]), 
which is independent of the overall normalization of f(a p ,t) (the pomeron flux). We see 
that there is a small (<5%) negative correction to the approximation fl22|) for a p and the 
same, positive, correction for y P . This effect can be understood intuitively as follows. 
The form of f(a p , t) favours low values of a p and therefore values of 90° < (p ep < 270°, i.e. 
cos0 ep < 0. In this region, the pomeron moves towards the virtual photon, thus increasing 

4 We use yfs = 296 GeV for all following numerical evaluations. 

5 We assume here that F 2 P is independent of t, i.e. that f(a p ,t) in (||) takes account of the full 
i-dependence. 
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the virtual photon energy yp 'seen' by the pomeron. Note that this effect decreases 
with increasing Q 2 and so will vanish in the asymptotic scaling limit. Furthermore, the 
deviation is proportional to the intercept 2a (0) — 1, which results in corrections of up to 



8% for 'hard' parametrizations of f(a p ,t), as suggested from BFKL phenomenology \'\ I 



In order to examine the factorizability of the diffractive structure function (^), we 
return to Eq. ([[J) in its fully differential form: 

da p dtdxdzdQ 2 ~ xQ A I 2 2 1 ' W ' p ' j 

-yl F E s (x,Q 2 ,a p ,t)} J^^5(z-z(x,Q 2 ,a p ,t)), (28) 
where we have made the replacement 

R x,Q ,a p ,t F 2 (x,Q 2 ,a p ,t) - F L (x,Q 2 ,a p ,t)' 29 

Assuming that the factorization of the structure functions F 2 and Ft gives a valid ap- 
proximation for the factorizability of the cross section, this can be expressed as 

da DS _ Ana 2 ( 1 + (1 - y) 2 



da p dtdxdzdQ 2 xQ 



Ff(z,Q 2 ,t) 



4 

F[(z, Q 2 , t) }f(a p , t) fj ^6(z - z(x, Q 2 , a p , t)). (30) 

After a simple integration over and x, restricted to the kinematically allowed values 
of the latter for fixed a p and z, we obtain finally an expression for the differential cross 
section similar to (|2|): 

da DS Aira 2 / Atz 2 \^ f y 2 1 » . „, 

'1-lZ+ on I dp/ \ F[{z,Q 2 ,t)f{a p ,t). 



da p dtdzdQ 2 zQ* \ Q 2 J { ° 2[1 + R p (z,Q 2 ,t)} 

(31) 

Assuming F 2 P and F£ to be independent of t, we can estimate of the magnitude of the 
Jacobian factor, 

r° ( Uz 2 \~^ 

,,,, „ 2 , /A 4 t ,n-\ ( o2 Igy^ 1+ gr ^ 

{J){z,Q ,a p ) = { 1 + — ^- )(z,Q,a p ) = , (32) 



Q 2 



dt f(a p ,t) 



which is shown in Fig. 0(b) as a function of z and Q 2 .[] We see that this Jacobian 
factor differs by less than 5% from unity for the whole kinematical range experimentally 

6 The a p dependence turns out to be negligible. 
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accessible at HERA. Together with the sytematic difference between yp and y, which is of 
about the same order, we find that the cross sections defined by (0) and (|3Tj) agree within 
a maximum deviation of 10%, which is attained only in the large- z region. For values of 
z < 0.4 the agreement is already better than 5%. Furthermore, both expressions become 
equal in the scaling limit Q 2 — > oo. As the experimental errors on the diffractive structure 
function are still well above these corrections ||14fl , and uncertainties arising from the R- 



factor are twice as big as these corrections, it seems appropriate at this time to factorize 
the diffractive structure function into a pomeron emission factor and a deep inelastic 
structure function of the pomeron, Eq. (H). When, in the future, the data improve and 
the full pomeron kinematics can be reconstructed, it should be kept in mind that this 
factorization is only an approximation to the factorization of the diffractive cross section, 
Eq. (§. 

A final point concerns the measured intercept of the pomeron trajectory. As a mea- 
surement of t is not possible at present, only an 'average' coupling of the pomeron to the 
proton can be determined: 

f(a p )= f dt f( a ,t) ~ a;" (eff) . (33) 

J — oo 

Using the DL-parametrization (|5|) for f(a p ), one finds n(eff) ~ 1.09 ±0.02. The error here 
represents the spread in n(eff) values as a p varies over the range 10~ 4 < a p < 10~ 2 . If 
a p is approximated by x/z, the effective power increases slightly to n(eff) ~ 1.11 ± 0.03, 
which is a non-negligible shift. Both these values are significantly lower than the 'naive' 
approximation n(eff) ~ 1 — 2a(0) = 1.17, and therefore this effect should be taken into 
account in comparing the measured intercept with model predictions. 

In summary, we have shown in this section that effects arising from an incomplete 
reconstruction of the pomeron kinematics at HERA give systematic corrections of only a 
few percent to a p , yp and the measured intercept of the pomeron trajectory. Furthermore 
we have demonstrated that the factorization of the diffractive structure function gives a 
correct approximation of the factorization of the diffractive cross section up to a relatively 
minor error, which vanishes in the large Q 2 scaling limit. 



3 Final-state electron-proton correlations 

It should be clear from the above discussion that identification of the scattered proton 
and measurement of its four momentum p' will provide a crucial test of the pomeron 
picture. In principle, this would allow a direct measurement of the parameters a p and t 
and hence of the pomeron emission factor /. However in practice it will be difficult to 
make a precision measurement of the proton energy, which would be needed to obtain 
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sufficient experimental resolution on a p and hence a precise determination of t. In the 
short term, it therefore seems more promising to test the t-dependence of / by using the 
angular correlation between the transverse momenta of the outgoing electron and proton, 
FigJ. 

As discussed in the previous section, the a p dependence of any Regge-motivated 
f(a p ,t) favours low values of a p , and therefore final state configurations in which the 
scattered electron and proton are approximately back-to-back. This effect will be en- 
hanced with increasing transverse momentum of the pomeron. Thus the distribution of 
events in the relative azimuthal angle <p ep is a measure of the average scale of t involved 
in the process. The <p ep dependence of the diffractive cross section can be parametrized 
in the form of a distribution function: 



x / Q2 v /3^ 

dt f I - + 2x\ 1 — cos <p ep , t 

dN „ Jtmin \z V xs Q i 

—(a;, z, Q\ <\> ep ) = * ro y^— x L (34) 



x 



where the lower limit on t arises from the physical range of the fractional proton momen- 
tum carried by the pomeron < a p < 1:[| 



Q 2 



90° < <p ep < 270° 



1 4z 2 (1 - Q 2 /xs) cos 2 <j) ep (r .r\ 

Imin — \ q2 /11\ 2 ^ ' 

— j— : K < 90°, ep > 270°. 



4(1 — Q 2 /xs) COS 2 (f) ep \x 

In practice, these bounds on t have minimal impact on the dN/d<p ep distribution, since 
one expects / to be strongly suppressed for \t\ values larger than a 'typical' hadronic scale 
ofC(lGeV 2 ). 

Figure § shows the predicted correlation between the outgoing electron and the rem- 
nant proton as a function of x, z and Q 2 . In fact it turns out that this function is almost 
independent of the ratio x/z, the naive expectation for a p . As expected from fl34|), the 
maximum asymmetry between the same-side and opposite-side hemispheres is obtained 
for low values of Q 2 and high values of z. Note that the effect reaches a magnitude of 
up to 30% for realistic HERA kinematics (Q 2 = 8GeV 2 ,z = 0.6), and hence should be 
distinguishable from statistical fluctuations. 

As we have discussed in detail in the previous section, the discrepancy between fac- 
torization at the level of the diffractive structure function and the diffractive cross section 
is of order —t/Q 2 , which is subleading to the y^i/Q dependence in ([34]). It is therefore 



7 This constraint is not to be confused with the more restrictive experimental cuts on the quantity 
x/z, since x and z are fixed in this distribution. 
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appropriate to use this angular distribution in connection with the factorized structure 
function (f|). Assuming the structure functions F[ and F 2 P to be independent of t, this 
yields the following result for the diffractive cross section: 



d ° W _ ^ U .,, V 2 ) f p (z Q*)±^L (x z O 2 



'ep 

(36) 

The error implicit in this expression due to the neglect of the Jacobian factor discussed 
in the previous section affects the normalization of dN/d<p ep , and leads to 



However this deviation is less that 5% for the kinematical range at HERA, since it only 
reparametrizes the Jacobian factor (^), which is small compared to the angular asym- 
metry of up to 30%. 

Eq. fl3"Bj ) can be used to extract the dN/ d<p ep distribution from the HERA data, since 
it only requires information on the coordinates of the remnant proton, and not on its 
momentum. This distribution can provide a crucial test of the applicability of DL-like 
parametrizations of f(a p , t). Furthermore, any t-dependence of F 2 P would result in devia- 
tions from the predicted z-dependence of dN/ d(p ep . In particular, a significant t-dependent 
contribution to F p would map the ^-dependence of F p onto the ^-dependence of dN/ d<p ep . 

4 Predictions for Ff and Ff 5 

4.1 Models for the partonic content of the pomeron 

The type and distribution of the parton constituents of the pomeron has been a topic 



of some debate Hl5[ . On one hand, it seems natural to assume that the pomeron is 



predominantly 'gluonic' []nj. On the other hand, the pomeron must couple to quarks at 
some level. In fact in Ref. || Donnachie and Landshoff have presented a prediction for 
the quark distribution in a pomeron 

zq p (z) « ±Cvrz(l - z), (38) 

with C ~ 0.17. This result is obtained from calculating the box diagram for 7*P — > qq, 
in the same way as the photon structure function is calculated in the parton model from 
the box diagram for 7*7 — > qq. A crucial difference for the above pomeron calculation 
is the softening of the pomeron-quark vertex by a form factor which suppresses large 
virtualities. This leads to the scaling behaviour (^) in the Q 2 — > 00 limit, in contrast to 
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the asymptotic growth g 7 (x, Q 2 ) ~ a(x) ln(<5 2 /A 2 ) obtained for the quark distributions 
in the photon. 

The absence of pointlike pomeron-quark couplings in the above model, which gives 
rise to asymptotic Bjorken scaling for the pomeron structure function, suggests that the 
partonic content of the pomeron is on a similar footing to that of any other hadron. In 
particular, we would expect the parton distributions to satisfy a momentum sum rule, 
fq + fg = 1, where f q (f g ) is the momentum fraction carried by quarks (gluons). If we 
take the Donnachie-Landshoff form (^) and assume three light flavours of quarks and 
antiquarks, we find f q = Cn/3 = 0.18 <C 1, which in turn suggests g p 3> q p . This is the 
basis of the model of Ingelman and Schlein || . If the gluon and quark distributions in the 
pomeron are both considered to be hard, as in Eq. for example, then good agreement 
is obtained with the shape of the UA8 jet rapidity distribution ||. However, preliminary 
analysis of the UA8 data has shown that the normalization of the data appears to be a 
factor of about 2 — 5 smaller (depending on the assumed mixture of quarks and gluons) 
than the prediction obtained from folding momentum-sum-rule constrained quark and 
gluon distributions with pomeron emission factors such as that given by Eqs. (|5]||) 0. 
There are two possible resolutions to this discrepancy. 

(i) It has been argued fllT|, [L8|, [L9| that there is an ambiguity in the definition of the 



pomeron emission factor f(a p ,t). Since the pomeron is not a 'real' particle, this 
factor cannot be interpreted as a probability in the strict sense, and therefore its nor- 
malization is process dependent. The normalization of the 'soft' / differs from the 
normalization of / in any hard process. Even the normalizations in different hard 
processes do not necessarily have to be the same. In the expression for the diffractive 
structure function (f|), one should therefore either replace the factor / (obtained by a 
certain procedure from 'soft' hadronic cross sections) by J\ff, or equivalently include 
an additional normalization factor J\f in the relationship between the pomeron struc- 
ture function and the quark distributions, i.e. F p = AfJ2 q e qQ F '■ We shall adopt the 
latter procedure below. 

(ii) The pomeron flux factor is (approximately) universal, but the pomeron has a point- 
like coupling to quarks, which gives rise to a direct contribution to the pomeron 
structure function. The analogy here is with the photon structure function, which 
grows asymptotically as q 1 ~ logQ 2 . A model of pomeron structure of this type, 
i.e. with q p ~ logQ 2 has been presented in Ref. |20"f . 



We would like to propose two very simple, physically motivated models for the pomeron's 
parton structure based on each of the above assumptions. It will turn out that at present 
both are able to give an adequate description of the diffractive structure function mea- 
sured by HI, but future more precise measurements, particularly of the Q 2 dependence, 
should be able to distinguish between them. 
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4.1.1 Model I 



We assume that at some bound-state scale, Ql = 2 GeV 2 (corresponding roughly to the 
mass scale of the glueball candidate reported in [Q), the pomeron is composed of valence 
gluons accompanied by a small amount of valence quarks and antiquarks.0 As the pomeron 
carries the quantum numbers of the vacuum, its quark and antiquark distributions have 
to be identical. Therefore, one only has to consider two parton distributions in the 
pomeron, the quark singlet S p = J2i(q[ + qC) an d the gluon. These are assumed to have 
the following, valence-like shapes at Q%: 



zH p (z, Ql) = f q {Ql) 6z(l - z), zg p (z, Q 2 ) = f g (Q 2 ) 6z(l - z) 



(39) 



For Q 2 > Ql additional quarks are generated dynamically, according to the GLAP 
evolution equations of QCD p2[ , and acquire a growing fraction of the pomeron's momen- 
tum. In fact, leading-order perturbative QCD predicts that the asymptotic momentum 
fractions are, regardless of the type of hadron, 



3n 



f 



16 + 3n 



16 



/ 



16 + 3n/ 



(40) 



Our model is also motivated by the success of the dynamical parton model for the proton 



structure functions [p3H , in which the proton is a mixture of valence-like quarks and 



gluons [24] at some low scale. 

In the evolution of these parton distributions we always define the quark singlet to 
be the sum of only three light quark flavours (u, d, s); contributions of heavy quarks to 
F p , of which we will only consider the dominant charm contribution, are incorporated 
by projecting the massive contribution from the 7*g — ► cc fusion process onto F p . This 
treatment of heavy quark contributions to deep inelastic structure functions has been 
shown 



to be more reliable than the construction of intrinsic heavy-quark distributions 
in the hadron, which then evolve like massless partons above a certain threshold. As 
argued in Ref. [25[j, quark mass effects clearly remain relevant even at energies above the 
HERA range, which calls into question a massless resummation of these contributions. 

For completeness, we will briefly outline the QCD treatment of the light and heavy 
quark contributions to the pomeron structure, although it is identical to the procedure 
in Refs. |25 , ^6| . The parton distributions at higher Q 2 are determined from the leading- 



i A similar model, but with a softer gluon distribution, has been discussed in 
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orderf] GLAP evolution equations 



d ( q f(z,Q 2 ) \ _ f 1 * I P„ f P* f ) ( gf (£,Q 2 ) 



(41) 



in which we keep the number of massless flavours fixed at n f = 3 in the splitting functions 
Pij, while the number of active flavours in the running of a s (Q 2 ) is determined by the 
Q 2 scale. This procedure results in continuous parton distributions and couplings at each 
flavour threshold, while A^ D is matched at each threshold. We use A^Q D (nf = 4) = 
200 MeV. 

Assuming that SU(3) flavour symmetry is already established at Qq, the contribution 
of the light quarks flavours to F p is just the singlet distribution times a charge factor: 

F 2 PMs \z,Q 2 ) = 2 -zi:(z,Q 2 ). (42) 
The massive charm contribution arising from photon-gluon fusion takes the form 

Ff> M V9 = 2 Zq ?-M± f^C (*-,§) S %,„% (43) 

271 Jaz y \V Q 2 J 

with the kinematical bound a = 1 + Am 2 c /Q 2 and the LO coefficient function 
C(C,r) = i[C 2 + (l-C) 2 + 4C(l-3C)r-8CV]lnl±| 

+ ^[-l+8C(l-C)-4C(l-CR (44) 

where 

f3 2 = 1 - i^. (45) 

It has been shown in [25] that a mass factorization scale of [i 2 = Am 2 for the gluon 
distribution in the above expression is the most appropriate choice with regard to the 
perturbative stability of the expression. We will use m c = 1.5 GeV in our numerical 
evaluations presented below. The complete prediction for F 2 P is therefore 

F 2 P (z, Q 2 ) = lzX p (z, Q 2 ) + F p ^\z, Q 2 , m 2 c ). (46) 



9 Although the full next-to-leading order technology is available, we do not consider it to be appropriate 
in this case. In the extraction of the diffactive structure function from the experimental data the 
diffractive i?-factor was set to zero, which can be only consistently accommodated in a leading-order 
parton distribution model. 
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The above treatment of the charm contribution takes proper account of the threshold 
behaviour which, as we will see in Section [5], makes a significant contribution to the Q 2 
dependence of the structure function. 

Finally, to take into account the ambiguity in the normalization of the pomeron flux 
factor discussed above, we multiply the structure function fl46|) by an overall normalization 
factor A/", 

F 2 p (z,Q 2 )^MF 2 p (z,Q 2 ). (47) 

We shall see that M ~ 2 gives a good representation of the HI data. Note that this factor 
cannot be directly compared to the 'discrepancy factor' obtained in the UA8 jet analysis 
M, since hard diffractive hadron-hadron processes do not general possess the factorization 



properties which would lead to universality of the pomeron emission factor [18, 27, 28 . 



4.1.2 Model II 

In the second approach discussed above, the pomeron structure function has two dis- 
tinct contributions, a 'resolved' component which consists of hadron-like quark and gluon 
constituents, as in Model I, and a 'direct' component coming from an assumed pointlike 
coupling of the pomeron to the quarks. This is in exact analogy to the photon structure 
function, and there is consequently no overall momentum sum rule. The Donnachie- 
Landshoff flux factor is used without modification. 

The resolved contribution is described in terms of parton distributions evolving from 
valence-like distributions for light quarks and gluons at a bound state scale of Qq = 
2GeV 2 : 

zZ p (z, Ql) = f q (Q 2 ) 6*(1 - *), V>, Ql) = fM) 6*(1 - z). (48) 
The contribution of the light resolved quarks to F 2 P is then: 

F P(M, S ), r e, ( ^ Q2) = Q2) _ (49) 

The resolved charm contribution arising from the fusion of a resolved gluon with the 
photon takes the same form as in Model I (|43|) : 

F P ^ es -(z,Q 2 ,m 2 c ) = if (c) (^,Q 2 ,m c 2 ) (50) 

The direct coupling of the pomeron to quarks with an unknown coupling c gives rise 
to a pointlike contribution to F p from pomeron-photon fusion.^] Here we take only 
the nonresummed contribution from the corresponding parton-level subprocesses, as the 

10 The treatment presented here is similar to the discussion in Ref. p0|. 
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resummed distributions tend to overestimate the magnitude of the distributions in the 
high- z region. The light-quark contribution to the structure function is then 



-,P(u,d,s), dir. 



•Q 2 ) = 2£ e 



1 8^2 



c 2 z 



z 2 + (l- z) 2 



111 



Q 2 (l-z) 2 



m 2 q z 



-z^(z,Q 2 ). 



(51) 

Here we assume the mass regulator for all light quark masses to be m q = 0.3 GeV. The 
full contribution from charm quarks above their production threshold 



Q 2 > Am 



2 

C l 



takes the form 



p P(c),dir./^ n 2 _2 



z,Q 2 ,m 2 c )=2q 2 c 



47T 5 



2 r ( Q 2 

c C [z, — 
\ mi 



The structure function of the pomeron is therefore 



F 2 P (z,Q 2 



7 ,P(u,d,s),Tes. 



z,Q')+Fl 



P(c), res. 



z,Q 2 ,m 2 c ) 



+F P(u, d ,s),^.^ Q 2 ) + F P(c)Ar. { ^ Q 2 ; m 2)_ 



(52) 



(53) 



(54) 



We shall see below that both models can be adjusted to give satisfactory agreement 
with the recent HI data, but give significantly different predictions for the Q 2 -dependence 
of F 2 P and for the charm content of diffractive scattering events. 



4.2 Q 2 evolution of F§ s 

The assumption that _F 2 is factorizable into an emission and a DIS part (|j) implies that 
the Q 2 dependence of F® s arises entirely from F 2 P . Furthermore, if the parton interpre- 
tation of F 2 P is valid, then this Q 2 dependence should be given by the standard GLAP 
evolution equations ( fill) of perturbative QCD. The observation of this Q 2 dependence is 
an important test of the whole approach, in particular of the factorizability of diffractive 
scattering cross sections. 

The Q 2 evolution of the resolved contribution to F 2 P is given directly by the GLAP 
equations (f41~p. For F® s we must fold the results with the pomeron flux factor /. In 
particular we can define 'diffractive' parton distributions in the proton by 

/ q DS (x,Q 2 ,t) \ /•! / q p (z,Q 2 ,t) \ 

{ g"(x, Q 2 , t))=L dzdap { g p (z, Q 2 , t) ) mz - x ' a ^ (55) 

where we have used flT£|), dropping the small corrections due to finite t and M 2 effects. 
Taking d/d In Q 2 of both sides, and using the fact that the resolved pomeron parton 
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distributions satisfy the GLAP equation, Eq. (f4lf) , gives 

dlnQ 2 \ g DS (x,Q 2 ,t) J 2tt Jo y ' p 1 w;j \ g p (y',Q 2 ,t) J 

xf(a P ,t) 8(z - x/a p ) S(z - rjy'), (56) 

where [P(^)] is the 2x2 matrix of splitting functions. Introducing 1 = / dy5(y — y'a p ) 
and integrating over y' and z gives 

a (q * (M \_ a ,m f^y lv{x/y)] ( (57) 



which is the usual GLAP equation, but now for the resolved diffractive parton distribu- 
tions. 

If the pomeron has a direct coupling to quarks as discussed in Model II, this will result 
in a direct contibution to the diffractive structure function 

F 2 DS > dil -(x,Q 2 ,t) = f'dzdaX^-^Q 2 ) f(a P ,t)5(z-x/a p ). (58) 
Jo 

In the scaling limit Q 2 — ■> oo, this term will dominate the Q 2 -dependence of F 1 ? 8 



p 



(59) 

where the sum runs over all flavours above their production threshold (|52[) . 

Therefore, one should find that the Q 2 dependence of both FJP s (x, Q 2 , t) andF^iz, Q 2 , t) 
is consistent with perturbative QCD while the corresponding parton distributions are re- 
lated by Eq. (|55|). A possible direct contribution to F 2 P will dominate the Q 2 -dependence 
of the diffractive structure function and should therefore be identifiable as the data im- 
prove. 

It is worth stressing that the Q 2 dependences of the proton structure function F 2 and 
F 2 at the same Bjorken x value are completely unrelated. In particular, F 2 rises rapidly 
with increasing Q 2 at small x as more and more slowly-moving partons are generated by 
the branching process. This rise is observed |2^] to be proportional to InQ 2 at fixed x, 
which is consistent with recent parametrizations [fj(| [30|, |31|] of the parton densities in the 



proton. In contrast, the resolved quarks in the pomeron are sampled at z values much 
larger than x, where the distributions evolve more slowly. 

For a pomeron without a direct coupling to quarks (Model I), one should therefore find 
that the fraction of diffractive events at fixed x is decreasing approximately like 1/ InQ 2 . 
A direct coupling (Model II) will result in a rise of F 2 DS proportional to InQ 2 at fixed 
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x. The intercept of this rise is unrelatedto the rise of the proton structure function and 
could in principle be used to determine the coupling strength of the pomeron to quarks 
from Eq. (|59|). 



5 Comparison with data 

The HI collaboration has recently measured [14] the diffractive structure function F 2 ^ 

d ° DS . . 4 ™ 2 i + (i-y) 2 F D(3 )( „ 2 v (m 

dxdQ 2 dx P - 2 F2 (AQ ' XP) (60) 

as a function of the three kinematic variables, {/3,Q 2 ,xp} where 

Q 2 







M 2 X + Q 



2 



x P = - « a P (61) 

with the approximations becoming exact when t = M 2 = 0. The variables oep, t and 
4> e p are not measured directly. By implication < <p ep < 2tt, and it is estimated that 



\t\ < 7 GeV 2 [0. 

By integrating Eq. ([!]) we obtain our prediction for the measured diffractive structure 
function 

F 2 D ^(P,Q 2 ,xp)= f'daJ^dtr^sLp---^) f(a p ,t)F 2 p (z,Q 2 ,t), 
Jo Jt min Jo 2n \ z Q l J 

(62) 

with z given in terms of the other variables by flT8|). Ignoring the t dependence every- 
where except in /, and setting the proton mass to zero, we obtain the simple factorizing 
approximation 



F 2 D{3 \(3,Q 2 ,x P ) 



max 



dt f{x P ,t) 



mm 



F 2 P ((3,Q 2 ), (63) 



which implies that the dependence of the structure function on ip should be universal, 
i.e. independent of (3 and Q 2 . Furthermore, if we substitute for / using Eq. fl5|) we find 



F 2 D ^((3,Q 2 ,x P ) « K xp n F 2 p (f3,Q 2 ). (64) 



Precisely this behaviour has recently been observed by the HI collaboration ||14j| . In fact 
their measured 'universal' power n of xp is n = 1.19 ± 0.06(stat.) ± 0. 07 (sys.), which 
is in excellent agreement with our prediction of 1.11 ± 0.03 (see Section 2) based on a 
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correct treatment of kinematics and using the pomeron emission factor of Donnachie and 
Landshoff §. 

The HI collaboration have also attempted to measure the pomeron structure function 
directly, by defining an xp-integrated diflractive structure function 

F 2 D (P,Q 2 ) = / dx P F? {3 \p,Q 2 ,x P ), (65) 

J0.0003 

where the range of integration is chosen to span the entire xp measurement ranges- 
According to the simple factorization hypothesis, the diflractive structure function is 
proportional to the pomeron structure function: 



with 



F 2 D (P,Q 2 )*AF 2 p (f3,Q 2 ), (66) 

i-0.05 fW 

A — da p dt f(a p ,t) « 1.5. (67) 

J0.0003 Jt min 

The numerical value in flB?|) corresponds to the DL form for /. In what follows we will use 



Eq. (|66D with A = 1.5 to convert the measured structure function |L4j into the pomeron 
structure function. 

In Ref. 0, data on F 2 D (/?, Q 2 ) are presented in four Q 2 bins, Q 2 = 8.5, 12, 25, 50 GeV 2 . 
In the first of these, the charm contribution should be relatively small, and hence F 2 D (/3, Q 2 ) 
can be directly compared with the predictions of our simple models for the light quark 
distributions, which allows us to tune the parameters of both models. 



5.1 Model I 

As the pomeron fulfils a momentum sum rule in this model, the first moment of F 2 , which 
is related to the the momentum fraction carried by quarks in the pomeron, contains the 
information needed to adjust the main parameters. Neglecting the charm contribution to 
F 2 P , this relation reads: 

A- 1 £ d(3 F 2 D ([3, Q 2 ) « £ d(3 F p ((3, Q 2 ) « 2 -N £ dz zE p (z, Q 2 ) = 2 -Mf q {Q 2 ). (68) 

The parameters M and f q (Ql) are strongly correlated — their product is essentially 
determined by the first moment of F 2 in the lowest Q 2 bin. Allowing for only a small 

11 Very recently a similar measurement has also been made by the ZEUS collaboration [3^| . The (3, Q 2 
dependence is consistent with the HI data. The normalization is different, however, because of the 
different x P range (0.00063 < x P < 0.01), which results in A(ZEUS)« 0.8. 



19 



charm contribution in this bin, we find the best agreement with the data for M = 2 and 



}2. 



the following momentum fractions of quarks and gluons at Qq. 

/ ff (Qg)=0.17, = 0.83. (69) 

Fig. [| shows the values of f q extracted from the HI data [14] in this way0 at the four Q 2 
values. Note that in the measured Q 2 range, the momentum fractions are predicted to 
vary only slightly with Q 2 . The apparent rise in the data has a simple interpretation as 
the onset of the charm contribution, as predicted by flUf) . 

In Fig. |5]the predictions of this model for the pomeron structure function are compared 
with the data, as defined by (pB|). The solid curves show the full prediction including the 
charm contribution, and the dotted curves are the contributions from the three light 
quarks only. We note that for this model, 

(i) an accurate description of the data without a direct pomeron coupling seems to be 
possible if the pomeron flux is adjusted; 

(ii) the variation of the dotted curves with Q 2 shows that the scaling violations predicted 
by the QCD evolution equations are rather weak in this kinematic range; 

(iii) the charm contribution grows rapidly above threshold (in fact, this growth is evi- 
dently responsible for the bulk of the predicted Q 2 dependence), and constitutes a 
significant fraction of the structure function at high Q 2 and low z; 

(iv) as Q 2 is increased to higher values, the pomeron structure function is expected to 
rise rapidly at low z and to decrease slowly at high z. 

Finally, in Fig. |6] we show the gluon and singlet (light) quark distributions in the 
pomeron, as predicted in this model.0 Since we are assuming exact SU(3) flavour sym- 
metry, the individual quark or antiquark distributions are simply q p = |S P . Note that as 
Q 2 increases, both the quark and gluon distributions evolve slowly to small z, as expected. 
The emergence of a small- ,2 'sea' of qq pairs can be seen at high Q 2 . 

12 The /3-integrated structure function in ( |68| ) is estimated by assuming that the structure function is 
independent of j3 at each Q 2 value. This is a very crude procedure, and we have no way of estimating 
the errors on the integral obtained by this method. Our comparison is therefore only semi-quantitative 
at best. 

13 The FORTRAN code for the distributions and structure functions in both models is available by elec- 
tronic mail from T.K.Gehrmann@durham.ac.uk 
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5.2 Model II 



In order to adjust the free parameters of this model (the ratio of resolved quarks and 
resolved gluons at Ql and the strength of the direct coupling c), we use the fact that the 
direct contribution (|5lD is concentrated in the high-z region. Therefore, the z = 0.065 
data points can be used to adjust the resolved contribution to the pomeron structure 
function. Best agreement is found for 

/ ff (Qg) = 0.3, = 0.7. (70) 

Having fixed the resolved contribution, the quark-pomeron coupling constant is deter- 
mined from the data to be c = 1. The parameters of this model are almost identical to 
the parameters of a similar model described in Ref . |2(| , obtained from an analysis of the 
rapidity and transverse energy distributions in diffractive photoproduction. 

Figure [7] shows the predictions of this model for the pomeron structure function com- 



pared with the HI data. The solid line represents the full F 2 , Eq. (|54f) , the dotted line 
is the sum of the resolved ( pEUp and direct fl5"T| ) light flavour contributions and the dashed 
line gives the resolved contribution summed over light ([49]) and charm (|50"D quarks. We 
note that for this model, 

(i) a good description of the HI data is obtained — in particular, the agreement in the 
Q 2 = 25 GeV 2 bin is better than for Model I; 

(ii) the pomeron structure function grows with Q 2 over the whole range in z, due to the 
dynamical generation of a qq sea at low z and the P — > qq splitting at high z; 

(iii) the charm contribution grows rapidly above threshold; in contrast to Model I, where 
the bulk of the charm events is concentrated at low z, we find a sizeable charm 
contribution in the high- z region. 

The different charm distribution in both models can be easily understood from the 
analytic forms of F% (Eqs. fl43|), fl50"D and (|53"D ): the coefficient function for the fusion of 
a virtual photon with a spin-1 boson to produce a massive quark-antiquark pair (53) is the 



same (up to colour and normalization factors) in all three contributions. In Model I and 
for the resolved contribution to Model II, this coefficient function is convoluted over the 
distribution of gluons in the pomeron. As this distribution vanishes for z — > 1, one will 
only find a small charm contribution directly below the threshold. The charm production 
from pomeron-photon fusion in Model II involves an (implicit) convolution of the above 
coefficient function with the distribution of 'pomerons in the pomeron', which is 5(1 — z) 
at lowest order. Due to the form of the coefficient function, this contribution is increasing 
almost up to the threshold, where it falls off steeply. 
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The ratio between the light quark and gluon densities is twice as big in Model II as it 
is in Model I. Therefore, we expect only a small fraction of resolved events in Model II to 
contain charm quarks, and the main charm contribution will come from the direct fusion 
process. The entirely different shapes of these contributions close to the threshold explain 
the difference between the charm distributions in the two models. 

Figure |8| shows the parton distributions in Model II, where again exact SU(3) flavour 
symmetry is assumed. The separate plots for the direct and resolved light-quark densities 
show the origin of the high- and low-z peaks in the quark singlet distribution. 

As both models give drastically different predictions for the Q 2 -evolution of the pomeron 
structure function in the high-z region and for the kinematical distribution of charm 
quarks, we expect that improved data will be able to distinguish between the different 
approaches. 

6 Conclusions 

The idea that the pomeron has partonic structure || has been given strong support by 
the recent measurements of the diffractive structure function at HERA. In this paper 
we have presented a detailed study of deep inelastic electron-pomeron scattering. We 
first derived the complete set of kinematic variables for the deep inelastic diffractive 
cross section. We showed that when expressed in terms of appropriate variables this 
cross section is expected to factorize into a pomeron structure function multiplied by a 
pomeron emission factor, the latter being obtainable from hadron-hadron cross sections.^ 
At present the variables which define the pomeron momentum are not directly measured, 
although they can be inferred from the observed hadronic final state. However, in terms of 
the measured variables the factorization is only approximate. In Section Q we quantified 
the corresponding systematic error, and showed that it was below the present level of 
experimental precision. 

When the remnant protons are eventually detected at HERA, it should be possi- 
ble to measure their scattering angle relative to the electron in the transverse plane. 
If the electron-pomeron scattering picture is correct, this distribution is predicted to 
be non-uniform, with a preference for back-to-back scattering. We presented quantita- 
tive predictions for this angular distribution in Sec. |3|, using the Donnachie-Landshoff 
parametrization for the pomeron emission factor. 

Finally, we presented two simple phenomenological models for the pomeron structure 
function. The first is based on the idea that at a 'bound-state' Q 2 scale, the pomeron con- 
sists predominantly of valence-like gluons, with a small admixture of valence-like quarks. 

14 Note that factorization is a 'high-Q 2 ' phenomenon, and will of course break down in the Q 2 — > 
limit, see for example [ ^8| . 
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A momentum sum rule is imposed. At higher Q 2 -scales the distributions are determined 
by standard GLAP perturbative evolution. Our starting quark distributions are identical 
in shape, and similar in size, to those calculated by Donnachie and Landshoff. In this 
model it is necessary to rescale the pomeron flux factor (by a factor of approximately 2 in 
the case of the DL function) to account for the normalization of the HI data, whereupon 
good agreement is obtained with the measured z and Q 2 dependence of the pomeron 
structure function. The light (u, d, s) quarks carry about 17-25% of the pomeron's mo- 
mentum in the range of Q 2 currently measured by HI. Note that the fact that our quark 
and gluon distributions are 'hard' in the HERA Q 2 range means that standard linear 
GLAP evolution should be perfectly adequate. Gluon-recombination effects, giving rise 
to non-linear evolution (as studied for example in Ref. ||), would eventually be expected 
to become important at very high Q 2 when the distributions have evolved to low z. 

Our second model allows for a point-like coupling of the pomeron to quarks, which 
generates an additional 'direct' component for the structure function, in analogy to the 
photon structure function. In this model the pomeron flux factor obtained from soft 
hadronic processes requires no further rescaling, since the direct component is large and 
positive. In contrast to the resolved (hadron-like) contribution, the direct contribution 
grows linearly with InQ 2 , at all z. Thus the asymptotic Q 2 —>■ oo behaviours of the two 
models are very different, even though the current data are not yet precise enough to 
discriminate between them. 

The experimentally-measured (z, Q 2 ) range of the pomeron structure function includes 
the charm quark threshold region. This requires special treatment, since the charm con- 
tribution to Fjf is expected to be significant above threshold. Motivated by the successful 
treatment of the charm content of the proton, we have calculated this using the photon- 
gluon fusion process, which takes the threshold kinematics correctly into account. We 
have found that the charm contribution to F% is indeed sizeable, especially at high Q 2 
and low z. Indeed, in Model I the rapid increase of the charm contribution with increasing 
Q 2 appears to account for the bulk of the observed Q 2 dependence. In Model II a similarly 
large charm contribution is obtained, but this time mainly from the direct coupling of the 
pomeron to charm quarks. The distinguishing feature between the two models is that in 
Model II the charm contirbution is more concentrated at large z (compare Figs. |5| and |7j). 

Our results on the quark and gluon content of the pomeron have many implications. 
As already mentioned, we expect that a significant fraction of hard diffractive scatter- 
ing events will contain charm, and our distributions provide a way of quantifying this. 
Especially in Model I, the overall magnitude of the gluon distribution compared to the 
quark distribution also predicts a large value for the pomeron's i?-factor. In particular, 
we expect R p ~ 0(1), in contrast to R ~ 0(a s ) <C 1 for the proton, which results in 
a similar magnitude of R DS . However, a consistent estimate of this would require a full 
next-to-leading order perturbative calculation, which is beyond the scope of the present 
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paper. 

In summary, we have shown that simple quark and gluon parton models of the 
pomeron, combined with pomeron emission factors extracted from soft hadronic pro- 
cesses, give an excellent description of the HI data. There are many ways in which this 
simple picture can be tested, and our two models distinguished, both at HERA and else- 
where. In the short term, the measurement of the <p ep correlation, the Q 2 dependence 
of the pomeron structure function particularly at large z, and the identification of the 
predicted large charm contribution to the diffractive structure function appear to offer 
the best possibilities. 
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Figure 1: Kinematics of deep inelastic electron-pomeron scattering 
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Figure 2: Systematic deviations after averaging over t and <f) ep , using the DL- 
parametrization for f(a p ,t) (|5|): (a) systematic relative deviation between a p and its 
approximation x/z as a function of x. The upper lines correspond to z = 0.2, the lower 
ones to z = 0.8. yp and y show the same systematic deviations with the opposite sign; 
(b) magnitude of the Jacobian factor defined in Eq. d32|) . 
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dN / 6<\> ep {x,z,Q ,(|) ) - distribution 
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Figure 3: dN/d<p ep (x, z, Q 2 , 4> ep ) distribution for fixed values of x/z and Q 2 , at z — 0.2 
(solid lines), z = 0.4 (dashed lines), z = 0.6 (dotted lines), and z = 0.8 (dot-dashed lines). 
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Fractions of total pomeron momentum (Model I) 
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Figure 4: Fractions of total pomeron momentum carried by light quarks and gluons 
as predicted by leading-order GLAP evolution for three light flavours in Model I (see 
text). The HI datapoints shown in comparison are the values for the momentum fraction 
carried by the sum of all light quarks under the naive assumption of a negligible direct 
charm contribution to F<[. A normalization of M = 2 for the pomeron emission factor of 
Donnachie and LandshofT is assumed. 
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Structure function of the pomeron (Model I) 
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Figure 5: The deep inelastic structure function of the pomeron F^iz, Q 2 ) constructed 
from the parton distributions in Model I (described in the text) of the three light-quark 
flavours (dotted line) and with an additional charm contribution from the photon-gluon 
fusion process (solid line). The HI data are obtained from values for the diffractive 
structure function in terms of these variables fL4| , divided by a pomeron emission factor 
of 1.5 (derived from the model of Donnachie and Landshoff). The theoretical predictions 
are scaled by a factor Af = 2, as discussed in the text. 
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Parton distributions in the pomeron (Model I) 



zg(z,Q 2 ) 



2GeV 2 
10GeV 2 
100 GeV 2 
1000 GeV 2 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



0.6 



0.4 



z£(z,Q 2 ) 



2 GeV 2 
10 GeV 2 
100 GeV 2 
1000 GeV 2 




Figure 6: Parton distributions in the pomeron, assuming a valence-like structure at Qq = 
2 GeV 2 and no pointlike coupling of the pomeron to quarks (Model I). The relative 
normalizations are chosen such that gluons carry 83% and light quarks carry 17% of the 
pomeron's momentum at Q%. 
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Structure function of the pomeron (Model II) 
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Figure 7: The deep inelastic structure function of the pomeron F^z.Q 2 ) constructed 
from the parton distributions in Model II (described in the text) for three light-quark 
flavours (dotted line) and with an additional charm contribution from the photon-gluon 
and photon-pomeron fusion processes (solid line). The dashed line repesents the resolved 
contribution to _F 2 P for all four quark flavours. The HI data are obtained from values for 
the diffractive structure function in terms of these variables |TJ]], divided by a pomeron 
emission factor of 1.5 (derived from the model of Donnachie and Landshoff). This emission 
factor is taken literally to be the pomeron flux factor. 
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Parton distributions in the pomeron (Model II) 
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Figure 8: Parton distributions in the pomeron, assuming a valence-like structure of the 
resolved component at = 2 GeV 2 and a pointlike coupling with strength c = 1 of the 
pomeron to quarks (Model II). The relative normalizations between resolved gluons and 
resolved light quarks are chosen to be 70:30 at Qq. The lower plots show the resolved 
(Eq. [49]) and direct (Eq. [51]) contributions to the light-quark singlet distribution. 
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